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ABSTRACT 

The accretion of hydrogen-rich matter onto C/0 and 0/Ne white dwarfs in binary systems leads 
to unstable thermonuclear ignition of the accreted envelope, triggering a convective thermonuclear 
runaway and a subsequent classical, recurrent, or symbiotic nova. Prompted by uncertainties in the 
composition at the base of the accreted envelope at the onset of convection, as well as the range of 
abundances detected in nova ejecta, we examine the effects of varying the composition of the accreted 
material. For high accretion rates and carbon mass fractions < 2 x 10"'^, we find that carbon, which is 
usually assumed to trigger the runaway via proton captures, is instead depleted and converted to ^*N. 
Additionally, we quantify the importance of •^He, finding that convection is triggered by ^He+^He 
reactions for '^He mass fractions > 2 x 10""^. These different triggering mechanisms, which occur for 
critical abundances relevant to many nova systems, alter the amount of mass that is accreted prior to 
a nova, causing the nova rate to depend on accreted composition. Upcoming deep optical surveys such 
as Pan-STARRS-1, Pan-STARRS-4, and the Large Synoptic Survey Telescope may allow us to detect 
the dependence of nova rates on accreted composition. Furthermore, the burning and depletion of 
^He with a mass fraction of 10"'^, which is lower than necessary for triggering convection, still has an 
observable effect, resulting in a pre-outburst brightening in disk quiescence to > Lq and an increase 
in effective temperature to 6.5 x 10^ K for a 1.0 Mq white dwarf accreting at 10~^ Mq yr"'^. 

Subject headings: accretion, accretion disks — binaries: close — instabilities — novae, cataclysmic 
variables — nuclear reactions, nucleosynthesis, abundances — white dwarfs 
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L INTRODUCTION 

White dwarfs (WDs) in cataclysmic variable (CV) and 
symbiotic binary systems accrete H-rich material from 
main sequence and red giant donors, respectively, with 
accre tion rates of M = 10~^^ — 10~^ Mq yr~^ (,Warner 
Il995| ). As the accreted envelope gains mass, compres- 
sion of the material at the base of the layer leads to a 
temperature increase, eventually triggering H-burning. 
Depending on M, the WD mass, and the composition of 
the accreted material, the eventual o utcome will either 
be steady therrnally s t able H-burning (lSienkiewicz|[l975l. 
pM; 'Paczvnski"1983'; 'Fuiimoto"1982bl: iShen fc BildstenI 
2007; Nomoto ct al. 2007), or a conv ective therm onuclear 
runaway (|Starrfield et al.lll972tlFuiimotolll982af) . observ- 
able as a classical, recurrent, or symbiotic nova when the 
convective zone nears the WD photosphere. (Hereafter, 
we will refer to the onset of convection and the ignition 
of the nova interchangeably.) 

It was first recogniz ed from numerical simulations (e.g., 
IStarrfield et al.iri972l ) that the mass fraction of CNO iso- 
topes, ^cnOj in the burning layers during a nova out- 
burst must be 10 — 50 times higher than the solar value 
in order to eject sufficient mass and produce a light- 
curve similar to those observed for fast novae. Abun- 
dance measurements confirm this metal enhancement in 
the ejecta of fast n ovae (for observational sum maries, see 
iGehrz et al.lll998l and lHachisu fc Katoll2006( l. This de- 
gree of enrichment cannot be explained by accretion from 
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an evolved donor or nucleosynthesis during the nova out- 
burst (e.g.. Imuran fc Livio 1986), thus some mechanism 
of core dredge-up and mixing with the envelope must be 
invoked. 

There is on-going debate over the effectiveness of pro- 
posed mixing mechanisms, which we detail in ^ In 
particular, it is unclear whether convection is initiated 
above or below the core-envelope interface in C-poor 
or C-rich material. With this in mind, our work ex- 
amines the effect of the accreted composition on the 
pre-ignition characteristics of nova systems, assuming no 
CNO enrichment prior to convection. We calculate ig- 
nition masses and pre-ignition luminosities for models 
with a large range of metallicities centered around solar 
{Z = 0.1 — 5.0 Zq), non-zero ^He mass fractions {X^ = 



IQ-^ - 0.005), high accretion rates (M = 10"'=' - 3 x 
10"'^ Mq yr"^), and a mass range of M = 0.6- 1.35 Mq. 
Note that the eventual outcome for the models with a 
combination of high accretion rate and low WD mass 
is steady and st able H-bu rning I Sienkiewicz 1975, 198C 
Paczvn ski 1983: Fu iimotolll982U ' IShen fc BildstenI 12007 



Nomoto ct al. 2007,), so we omit these models from our 



results. 

In 12] and 21 we use analytic approximations to mo- 
tivate the more exact numerical study described in Sj5j 
Below a critical accreted carbon mass fraction X12 ~ 
2 X 10^"^, which is coincidentally near the solar value, 
we find that ^^C is depleted and converted to ^''N prior 
to unstable ignition, so that convection is triggered by 
proton captures onto ^''N. We also find that convection 
is triggered by ^He-|-'^He reactions for '^He mass fractions 
X3 > 2 X 10"'^. These different nova triggers change 
the pre-outburst luminosity of the nova system, which 
we detail in ^5.5^ The ignition mass is also affected by 
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the triggering mechanism, which results in a previously 
unconsidered dependence of galactic nova rates on com- 
position. This effect is especially relevant given the up- 
coming flood of data from optical transient surveys such 
as Pan-STARRS-1, Pan-STARRS-4, and the Large Syn- 
optic Survey Telescope, which will measure nova rates in 
external galaxies with greater accuracy than available in 
current data. We speculate on the observational conse- 
quences and summarize our work in ^ 

2. MOTIVATION FOR OUR STUDY AND 
JUSTIFICATION OF ASSUMPTIONS 

Core-envelope mixing models differ on the mech- 
anism by which core material is brought into the 
accre te d layer. Chemical di ffusio n (Prialnik & Kovctzf 
1984| iKovetz fc PrialnikI Il985t Prialnik & Kovetz 
T995I: lYaron et al.l l200fiD and shear mixing caused 
by differential rotation of the accreted material 



(iKippenhahn fc Thomai |1978t iLivio fc TruranI 119871 : 
iKutter &: Sparks! 1X989^ Ulexakis et al.l I2004D result 
in pre-convective penetration of a small amount of 
hydrogen into the underlying material and vice versa. If 
the material below the layer is C-rich, p+^^C reactions 
trigger convection, which homogenizes the envelope 
and the entrained core material. These pre-convective 
enrichment models differ from mechanisms in which 
convection is triggered above the core in accreted 
material whose composition is determined by that of 
the donor star and is thus relatively C-poor. The 
introduction of core material into the envelope for 
these convective enrichment models is caused by the 
convective m otion itself, eit her via convective overshoot 
into the core (jWooslevI 19861) or shear mixing induced by 
the c onvective eddies (jClasner et al.|[l997l iRosner et al.l 
[2OOI . 

As of yet, no mixing mechanism has definitively proven 
itself successful in explaining the enrichments of all 
novae . For example, the mu lti-c ycle diffusion stud- 
ies of lPrialnik fc Kovetzl ()1995| ) and lYaron eFan (120051 ) 
are initiated with matter accreted directly onto naked 
C/0 cores, yet the accreting WD m ay be 0/Ne in as 
many as 1/3 of all observed n ovae (Truran & Livio 1986; 
Ritter et al] 119911 : iLivio fc Truran 1994, : Gil-Pons et al 



20031) ■ For these systems, diffusion would not lead to 



the initiation of convection below the accreted layer and 
subsequent core dredge-up becau se the underlyiiig mate - 
rial is not C-rich. The studi es of ILivio Jg^Truran' l igSTl ). 
iFuiimotol (|1988l I1993D . and iPiro fc Bilds'ten (2004| ) rule 
out any significant differential rotation between the ac- 
creted layer and the core, which cast s doubt on acc r etion- 
induc ed shear mixing mechanisms. iKercek et al.l (|1998l . 
II999') find that convective overshoot and shearing do not 
sufficiently enrich the envelope to produce a fast nova, 
although possible probl ems with their boun dary condi- 
tions are pointed out bv lGlasner et al.l ()2005l ). Moreover, 
some recurrent nova, which are novae with recurrence 
times ~ 30 yr, do not show over- abundances of metals in 
their electa dWilliams et al.lll981 : WiUiams 1982 : Warner 



Il995t iHachisu fc Katoll2001l )7 possibly due to a large he- 
lium buffer above the core. 

Thus, it is unclear how much the envelope will be 
enriched in metals prior to nova ignition for M > 
10^^ Mq yr~^. Some nova studies that include the ac- 
cretion phase assume that the accreting envelope is pre- 



enriched by the co re and consists of up to 50% core mate- 
rial b y mass (e.g. JJose fc Hernandll998l : IStarrfield et al.l 
Il998f) . However, the uncertainties involved in the mix- 
ing mechanisms coupled with the lack of observed metal 
enrichment in recurrent novae lead us to examine the 
consequences of assuming no C-enrichment in the ac- 
creted envelope prior to the onset of convection. There 
are several pre vious studies that also foll ow this treat- 
ment (Starrfic ld et all I1985I. Il98 8. 2000 jTruran et all 
Il988t jPiersanti et al.l 120001 : iJose et all l2(507l ). but none 
sufficiently samples the full parame t er spa c e in which 
we are inter e sted: IStarrfield et"an (|1985L 11988) and 
iTruran et al.l ()1988l ) consider Z = 0.02 accretio i i onto 
M > 1.35 Mq WDs; the models of lStarrfield eTall (I2OOOI ) 
have metallicity Z = 1Q~^ or 0.02, as motivated by novae 
in the Large Magellanic Cloud, with lowe r accretion rates 
5. 10~ ^ Mq yr~^ than what we study; iPiersanti et al.l 
(|2000f ) examine accretion with 3 metallicities (Z = 0.02, 
10-3, and lO"'^) onto WDs wit h masses < ( 3 . 68 , lower 
than our parameter range; andl Jose et al.l (|2007f ). in their 
study of novae in primordial binaries, have only solar and 
very sub-solar metallicity 1.35 Mq models with Z — 0.02, 
2 X 10-6, or 10-'^, with Af = 2 x 10-^*^ Mq yr-^, which 
is lower than our range. Moreover, none of these studies 
consider the effect of '^He, which ca n play a dominant role 
in tri ggering the nova (|Sharaa980l:lTownslev fc BildstenI 
l2f)f)l . 

In our study, we assume that the effect of chemical 
diffusion is negligible. This assumption, and thus our 
results, are invalid if the material directly below the ac- 
creted layer is C-rich. However, as we have described 
above, many nova systems exist in which the underlying 
material is C-poor and diffusion is indeed negligible. Our 
results only apply to these systems. 

3. CARBON DEPLETION PRIOR TO UNSTABLE 
IGNITION 

If the accretion rate in a CV is lower than the minimum 
rate for stability, the result will be a hydrogen shell flash. 
In this section, we calculate the ignition conditions for 
these thermonuclear novae. Throughout this study, we 
make the assumption that the accreted layer is thin, with 
a pressure scale height at the envelope base, h = Ph/ p^g, 
much less than the WD core radius, R. The subscript b 
refers to the base of the accreted envelope, g = GM/ is 
the gravitational acceleration, assumed constant because 
h R, and M is the WD core mass. The ratio of the 
scale height to the WD radius for an ideal gas equation 
of state is the ratio of the thermal energy to gravitational 
energy. 



h 
R 



kBTb/p.mp 
GM/R 



Ml 



(1) 



where T7 is the base temperature in units of 10^ K, Mi 
is the WD core mass in units of Mq, Rg = R/10^ cm, 
the proton mass is nip, and the atomic mass per particle 
is fj, = 0.6 for solar composition. For typical ignition 
conditions, Tf, ~ 2 x 10^ K, so the shell is very thin for the 
entire accretion phase. The base pressure is independent 
of the temperature in the thin-shell limit and is 



Ph 



GMM^r 
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(2) 
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where M^nv is the envelope mass. This assumption wiU 
not be vahd once the temperature rises during the ther- 
monuclear runaway. For a 0.6 Mq WD, h ~ R when 
Tb ~ 7 X l{f K. 

We first estima te the lum inos ity in the accreting 
layer fo llowing Nomotd ()1982| ) and lTownslev~fc BildstenI 
(|2004[ ). When material accretes onto the WD surface, 
the gravitational energy, GM/R, is re leased and radi- 
ated by the spreading iDOundary layer (jPiro fc BildstenI 
[200l and is not carried into the star, because the ther- 
mal timescale at the photosphere for luminosities of order 
the accretion luminosity is far shorter than the accretion 
timescale. Instead, prior to the onset of nuclear burn- 
ing, the pre-ignition luminosity exiting the deep accret- 
ing layer is produced by entropy released as the material 
accumulates. The entropy equation yields the compres- 
sional luminosity at the surface, 



ds 

Ml T-dP, 



(3) 



where s is the specific entropy, and we have neglected 
the term ds/dt\p. The lower bound, P' , is the depth at 
which the thermal time is equal to the time for which 
accretion has been on-going, so that the luminosity pro- 
duced there has had time to make its way through the 
envelope. For illustration, we assume Kramers' opac- 
ity (k oc pT"^/^), ideal gas (P oc pT), and a con- 
stant luminosity above P', so that P{r)^ oc T(r)^^/^. 
An ideal gas has s = /cb In (T^^^/p) / finip, which yields 
ds/dP = —7kB/l7fimpP. This gives 



-t'comp — T-^^" 

4 fiirip 



0.4 Lp, M_ 



/ J., 



(4) 



where T' is the temperature at P', M_8 is the mass ac- 
cretion rate in units of 10~^ Mq yr~^, and we have set 
/Li = 0.6. If the opacity is due to electron scattering, the 
pre-f actor becomes 3/2 instead of 7/4, so the exact re- 
lation is only weakly dependent on the form of radiative 
opacity. 

The thermal time at P' is ithorm = cpT'Mo„v/^comp, 
where cp = bkn /'^H'm.p is the specific heat at con- 
stant pressure for an ideal gas, Mg^^ is the mass in the 
layer above P', and we use a one-zone app roximation, 
dL/dM - Lcomp/M' (e.g.. lPaczvn"skilll983h . The time 
to accrete an envelope mass M^nv is iacc = M^nv/M. To 
find the depth from which luminosity is able to escape 
during accretion, we set the two timescales equal and use 
equation yielding P' ~ P(,; i.e., most of the luminos- 
ity in the envelope comes from only the envelope itself, 
and so we neglect the com pressional luminosity from the 
core (see the appendix of iTownslev &: BildstenI [200l for 
further discussion of the core's role). Thus, the compres- 
sional luminosity is given by equation ([4]), with T' = Tb. 

Using the radiative diffusion equation with Kramers' 
opacity, 



P 



g cm 



-3 



-7/2 



(5) 



whe re — 10^^ cm^ g~^ from fitting to OPAL opaci- 
ties (jlglesias fc Rogerslll993l 11996.) for solar composition 
around T = 10^ K and p = 10^ g cm""^, we derive the 



temperature at the base of the accreting layer as a func 
tion of pb, 

/ • \ 2/11 



Tb = 1.37 X 10^ K 



(6) 



where p^ = pb/W"^ g cm"-^. We have assumed solar 
metallicity, but this result is nearly independent of com- 
position. The bottom of the layer follows the trajec- 
tory given by equation ([6]) until nuclear burning becomes 
comparable to compressional heating, i.e., when L„uc ~ 
^compi where Lnuc is the luminosity produced by nuclear 
burning. For high accretion rates > 10~^ Mq yr~^, H- 
burning occurs via CNO reactions when base conditions 
reach ~ 2 x 10^ K and pb ~ lO'^ g cm"*^ (ignoring 
■^He-burning). If the accreting material has near-solar 
isotopic ratios, the most relevant isotope is ^^C, since 
proton captures onto ^^N are slower than onto ^^C, and 
^^O docs not participate in the CNO cycle at these tem- 
peratures. Moreover, p + p reactions are unimportant 
at Tft ~ 2 X 10^ K because the lifetime of a proton 
with respect to self-burning is ~ 10 times longer than 
with respect to consumption by ^^C nuclei. Thus, pro- 
ton captures onto ^^C will be the first non- negligible re- 
action. These are quickly followed by the /3-decay of 
^■^N (with a half- life of Ti/2 = 603 s) and proton cap- 
tures onto ^^C (~ 4 times more rapid than onto ^^C), so 
that we approximate the first nuclear reactions of interest 
as the conversion of ^^C to ^"'N at the ^^C proton cap- 
ture rate. This reaction chain releases a specific energy 
X12E12 = 8.8 X 10i4(Xi2/1 0-^) erg g-^. 

Linear stability analysis (|Fuiimoto et al.l [19811 ) shows 
that nuclear burning is unstable in a constant-pressure 
thin shell if 



dT 



> 



decooi 



dT 



(7) 



or e 



nucXnuc 



> CcooiXcooh where Cnuc is the nuclear en- 
ergy generation rate, the one-zone approximation to the 
cooling rate is ecooi L/Mcnv, and x = dlne/dlnT\p. 
For Kramers' opacity and ideal gas, Xcooi = 17/2. The 
cooling rate is rewritten as 



^cool 



TkeTb 1 

4 jlTflp tg^Q(^ 



(8) 



after substituting the expression for compressional lumi- 
nosity (eq. [3]). 

It is often stated that the energy released by the 
conversion of ^■^C to "'^''N triggers the thermonuclear 
runaway for high M no vae ( e.g.. iFuiimotol Il982bl: 
Townslev fc BildstenI I200I [Josl I2Q05I: iJose fc Hernand 



2007f ). However, we show here that under some condi- 
tions this reaction will not release enough heat to trigger 
unstable ignition conditions before the ^^C is depleted. 
In this case, all available ^^C converts to ^"'N, which will 
ignite later at a larger pressure and temperature. Car- 
bon depletion occurs when the accretion and burning 
timescales are comparable, with the C-burning timescale 
given by ti2 = A"i2-Ei2/ei2, where £12 is the rate of energy 
generation from conversion of ^^C to ^^N. The condition 
for stable ^^C depletion is then 



^i2<3xl0- % ^i^. 

2 X 10" K 8.5 X12 



(9) 
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Fig. 1. — Analytic ^^C ignition and depletion curves for a 
1.0 A/q WD. Dotted lines demarcate the boundary between stable 
'p-^-^'^G burning (lower-left region) and unstable ignition {upper- 
right region) for the given accreted metallicities. The short-dashed 
line shows where tacc = ti2, which is independent of M and Z; 
above this line, ^■^C has been depleted. The solid lines follow 
the trajectories at the base of envelopes accreting at M = 10~^ 
and 10~* Mq yr"'^; as more material accumulates, the base be- 
comes hotter and denser, and the lines move to the upper-right, 
where carbon either ignites unstably or depletes. If carbon is de- 
pleted, the base continues to heat and compress until ^''N begins 
to burn and the full CNO is ignited (long-dashed line). To the left 
of the dashed-dotted line, the base is mildly degenerate or non- 
degenerate. 



Thus, for mass fractions below a critical value, coinciden- 
tally near the solar mass fraction of 2.2 x 10~^, carbon 
will deplete before triggering a nova. This value is cer- 
tainly relevant for novae in systems with low-metallicity 
donors. Furthermore, low ^^C/^'*N ratios can occur 
when mass transfer h as revealed a CNO-processed core 
(jSchenker et al.l [20021 and references therein). In these 
cases, the carbon mass fraction of the accreted material 
will be well below the solar value because proton captures 
onto ^^N are the slow step of the CNO cycle, and thus 
the donor's CNO nuclei are mostly in the form of ^''N, re- 
sulting in accreted carbon mass fractions ~ 10~^ or lower 
(jSchenker et al.ll2002f l. For high M CV systems, evolved 
donors such as these a re likely common. A population 
synthesis calculation bv lPodsiadlowski et al.l ([2003) finds 
that the majority of CVs with orbital periods Porb > 5 
h have an evolved secondary. Observationally, a study of 
UV line fl ux ratios of CVs both above and below the pe- 
riod gap (Gan sicke et al.ir2003f ) concludes that as much 
as 10 — 15% of their sample might have evolved donors 
with anomalously low ^^C/^^N abundance ratios due to 
CNO processing. 

Figure [1] shows conditions for carbon depletion {short- 
dashed line) and stable vs. unstable C-burning with 
varying accreted metallicities {dotted lines) for a 1.0 Mq 
WD. The region to the left of the dashed-dotted line has 
Fermi energy Ep < SkeT, so that it is mildly degen- 
erate or non-degenerate. Free-free opacity dominates in 
the region of interest. Also shown are the rising tem- 
perature and density at the base of accreting layers with 
M ~ 10^^ and 10~^ M© yr^^ {solid lines). For mass 



fractions lower than near-solar, carbon depletes before it 
ignites unstably and is instead converted to ^"^N. In these 
cases, the base of the layer continues to become hotter 
and denser until it reaches burning conditions for proton 
captures onto "'^'^N {long-dashed line). At the tempera- 
tures and densities corresponding to these M's, the CNO 
cycle as a whole is thermally unstable, and ignition is in- 
evitable. However, because the layer must be hotter and 
denser to burn the ^"^N, more material must be accreted 
before ignition than in the case of carbon ignition. 

4. THE SIGNIFICANCE OF ^HE 

The possible importa nce of "^He for novae was first 
studied by ,Schatzmannl (|195lD in the context of a the- 
ory of novae pow ered by thermonuclea r detonations. 
iMacDonaldl (|1983f l and llben fc Tutukwl (|1984f ) noted 
that the presence of '^He can decr ease t h e env elope 
mass nee ded for nov a ignition, and ISharal ()1980( ) and 
[ib w nsley fc BildstenI ([2004) looked more closely at its 
role in triggering the convective shell burning phase of a 
nova. Although these studies found that "^He can play a 
large role in the onset of a nova, its effects have not been 
quantified in detail. 

For typical nova conditions, ^He-burning takes place 
via ^He(^He, 2p)'*He. The next-fastest reaction that con- 
sumes '^He is ^He-|-*He, which has an unscreened reac- 
tion rate that is a factor of ~ 10"^ slower for X3 = 10~^, 
and thus ^He-|-^He is the only reaction to consider. To 
gauge the importance of '^He-burning for nova ignitions, 
we compare its relevant characteristics to proton cap- 
tures onto ^^C, the assumed trigger for nova systems with 
M > 10^^ Mq yr~^. The energy per "^He self-reaction is 
12.9 MeV, slightly larger than the 11.0 MeV released by 
the reaction chain that converts ^^C to "'^''N. Moreover, 
the '^He unscreened reaction timescale, i.e., the e-folding 
lifetime for ^He nuclei, defined for a generic nucleus i as 
rii/ti = drii/dti, is much shorter than that of ^^C. The 
ratio of timescales is 



tl2 Tr,3 ^3 / 6.56 

lO-" — exp 

t3 Ah \ Tr 



1/3 



(10) 



which is roughly a factor of 200 for 2 x 10^ K, A3 10"^^ 
and hydrogen mass fraction Ah = 0.75. Thus, for A3 > 
10~^, '^He nuclei will begin burning via self-reactions be- 
fore ^^C nuclei have a chance to capture protons. 

Since the '^He reaction releases a similar amount of en- 
ergy and has a similar temperature dependence to that 
of the ^^C proton capture, we must also consider the 
possibility of '^He depletion prior to nova ignition. An 
analysis like that of the previous section yields the same 
critical mass fraction ~ 3 x 10"'^ as in equation Again 
coincidentally, this critical mass fraction for •^He deple- 
tion is in the neighborhood of the expected value for 
mass-transferring binary systems. As mass loss uncovers 
deeper parts of the donor star, material that has been 
proc essed by H-burning can make its way to the sur- 
face. iD'Antona fc Mazzitellil (|1982f l and llben fc Tutukovl 
dl984l ) find that the mass fraction of accreted '^He can 
reach values as high as 4 x 10~^ during the evolution 
of systems with low-mass donors, and thus the possi- 
bility exists for ■^He to tri gger a nova prior to de plet- 
ing. For this reason, Townslev fc BildstenI (|2004[ ) in- 
cluded calculations with A3 = 0.001 and 0.005. For 
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TABLE 1 
Model Core Temperatures 



Mass [Mq] 


M [Mq yr-1] 


Tc [K] 


0.6 


10-s 


8.00 X 10^^ 


0.6 


io-» 


1.04 X 10^ 


0.6 


10-^ 


1.34 X 10^ 


0.6 


3 X 10"'^ 


1.51 X 10'' 


1.0 


io-« 


8.99 X 10'^ 


1.0 


io-» 


1.10 X 10^ 


1.0 


10-^ 


1.45 X 10^ 


1.0 


3 X lO"'^ 


1.70 X 10^ 


1.2 


10-9 


9.18 X 10^ 


1.2 


io-» 


1.20 X lO'^ 


1.2 


10-' 


1.80 X 10' 


1.2 


3 X 10"'' 


2.29 X 10' 


1.35 


10-9 


9.25 X 10« 


1.35 


10-s 


1.29 X lO'^ 


1.35 


10-^ 


2.18 X lO'^ 


1.35 


3 X lO-'^ 


3.00 X 10'^ 



M > 3 X 10"^" Mq yr~^, they found that ignitions were 
triggered by ^He, with the ignition mass Mjgn roughly de- 
creasing by a factor of 2 when increasing from 0.001 
to 0.005. 

5. NUMERICAL SIMULATION 

The approximations that we have made in the ana- 
lytic work of and 2] and the proximity of the criti- 
cal ^^C and ^He abundances to relevant solar and CV 
abundances motivate a more exact analysis. In partic- 
ular, electron degeneracy pressure, other opacities, the 
exchange of heat between the envelope and core, and the 
term ds/dt\p that is neglected in equation jS]) must be 
included. Moreover, the one-zone approximation of set- 
ting dL/dM ~ L/Menv is problematic, because burning 
occurs in a narrow layer and is certainly not a linear 
function of the whole atmosphere. 

5.1. Model 

We developed a time-dependent explicit Runge-Kutta 
code for a one-dimensional grid of 100 zones covering a 
pressure range of 3 x 10^^ — 3 x 10^° dyne cm"^. Since 
the base of the envelope at ignition is typically between 
3 X 10^^ and 3 x 10^^ dyne cm"^, this choice of zon- 
ing gives greater than a factor of 10 range in pressure 
above and below the region of interest. For the 1.35 Mq 
model with M = 10~^ Mq yr"^, the pressure bound- 
aries were changed to 10^^ — 10^^ dyne cm"^ because 
the accreted layer reached a depth that was too close to 
3 X 10^" dyne cm"^ . The zones are spaced logarithmically 
in pressure (thus there are roughly 33 zones per decade 
of pressure) to better resolve the accreted layer and to 
avoid over-resolution of the bottom zones. The layer is 
spherically symmetric, appropriate for the depth s of in- 
teres t if there is negligible differen tial rotation (Fu iimotol 
119881 Il993t IPiro fc BildstenI 12004 ). and plane-parallel, a 
good approximation as discussed in 521 

The thermal evolution o f the core during CN cycle s 
was considered in detail bv'Tow nsley fc BildstenI (|2004j ). 
who found that heating and cooling during the nova cycle 
results in an equilibrium core temperature, Tc^cq-* Thus, 

* CVs spend only a short amount of time at M > 10"* Mq yr" 



TABLE 2 

WD Core Properties for Tc = lO' K and Outer Boundary 
Condition P;, = 10^® dyne cm"^ 



Mass [Mq] 


Radius [cm] 


g [cm s-2] 


0.6 


8.39 X 10** 


1.13 X 10* 


1.0 


5.40 X 10* 


4.55 X 10* 


1.2 


3.83 X 10* 


1.09 X 10'-* 


1.35 


2.26 X 10* 


3.52 X 10^ 



for our models, the initial thermal profile prior to accre- 
tion is assumed to be a radiative-zero solution that gives 
a core temperature Tc = Tc. oq. The equilibrium Tc de- 
pends on the accreted composition, but, in order to limit 
the parameters of this study, a single representative Tc 
is used for each M and M n iodel; these are extrapolated 
from iTownslev fc BildstenI (|2004f ) and Townsley (2007, 
priv. comm.) and are shown in Table [1] Fortunately, 
the properties of the accreting layer are relatively insen- 
sitive to Tc for these high M's. For example, increasing 
Tc by a factor of 2 decreases Mjgn by only < 10%. 

At each time-step, the temperature for the top-most 
zone is set with respect to that of the zone directly below 
it according to a power-law solution obtained by assum- 
ing a radiative-zero atmosphere above our grid.^ The 
thermal boundary condition for the bottom-most zone 
is such that the flux there is equal to that of the zone 
directly above it. While locally incorrect, this bottom 
boundary condition has no effect on the thermal prop- 
erties of the accreted envelope as the thermal time at 
the bottom of our grid is much longer than the accre- 
tion time prior to ignition. Thus, any inaccuracies in the 
bottom-most zones will have no effect on the region of 
interest. 

The WD core structure is assumed to be constant dur- 
ing the nova cycle. This is an excellent approximation, 
as the WD's central pressure, 10^'^ — 10^^ dyne cm"^ for 
M = 0.6 — 1.35 Mq, is much larger than that of the 
accreted layer. The WD core radius is calculated for an 
isothermal core that is half ^^C and half ^^O by mass with 
outer boundary condition = 10^^ dyne cm"^. The ra- 
dius is relatively insensitive to Tc for Tc < 2 x 10*" K, so 
all radii are calculated with Tc = 10*" K. The resulting 
radii and gravitational accelerations for our 4 models are 
shown in Table H Calculations for 0/Ne WDs give the 
same radh to within ~ 1%, because ^^C, ^^O, and ^°Ne 
have the same charge-to-mass ratiO; 

The equation of state (ISaumon et alJ Il995l : 
iTimmes fc Swestvl l2000t IR ogers fc Navfonovl [200l) . 
opacitv (jlgl esias fc Rog er j[l993t 1 19961 ), nucle ar burning 
netw ork jTimmcs 1999), neutri no cooling (jltoh et al.l 
Il996f l. and electron screening (jOraboske et al.l 119731 : 

l|Howell et aLI20011) . and so it is unclear if the core temperature will 
have time to reach equilibrium before the system has evolved ap- 
preciably llEpelstain et a l.'2007i. However, as mentioned later, the 
results of our study are largely unaffected by the core temperature. 

^ Changing the upper boundary condition to a constant tem- 
perature has a negligible effect on the long-term evolution of the 
layer. This is unsurprising because, for the pressure at the top of 
the grid, Ptop, the ratio of the thermal time to the accretion time 
at the base is Ptop/Pfc. For the grid we use and typical ignition 
pressures, this ratio is < 0.1, and so any differences in thermal 
conditions at the top of the grid are radiated away and not carried 
deeper into the star. 
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TABLE 3 

Solar Photosphere Mass Fractions from ILoddersI (120031') 



Element 


Mass fraction 


iH 


0.749 


4He 


0.237 


12q 


2.21 X 10^^ 


14n 


0.71 X 10^^ 


16o 


5.87 X 10^^ 


*Nc 


1.12 X 10"^ 


24 Mg 


0.64 X 10"^ 



lAlastuev fc Jancovicil 119781 : lltoh et al.l Il979f) are cal- 
culated using the MESA code package.*^ The MESA 
basic nuclear network, which tracks the abundances 

of IH, "He, l^C, 14^^ 16q^ 20p^g^ g^^^ 24]y[g^ ^^^^ 

modified to explicitly follow the '^He reactions '^He+'^He 
and ^He+^He. The nuclear reactions consist of the p-p 
chains (including the pep reaction, although for our 
high M's > 1 0~^ M(7i yr"-'^, it typic ally has a negligible 
contribution; IStarrfield erall l2007t ). the CNO cycles, 
a-burning up to ^''Mg, and C/0 burning. Isotopes 
that are not explicitly tracked are assumed to have 
equilibrium abundances. 

Solar composition of the accreted material is defined 
by the recomm ended solar photosphere abundances in 
iLoddersI (|2003D . The relevant elemental abundances by 
mass fraction are shown in Table [3l note the reduction 
in solar met allicity as compared to previous studies, such 
as lGrevesse &: Sauv al (1998). For models with non-solar 
metallicities, the mass fractions of ^^C, ^'^'N, ^^O, ^°Ne, 
and ^^Mg are adjusted by the same multiplicative factor^ 
and the difference made up in Xh- Likewise, any non- 
zero '^He mass fraction is subtracted from Xh- 

A monotonic transport firs t -orde r advection scheme 
(|van Leed fl97l iHawlev et aT]ll984f ). modified for log- 
arithmic coordinates, is utilized to simulate accretion. 
While certainly more accurate than zeroth-order donor 
cell advection, it is still subject to numerical advection 
that smoothes out what should be a step-function ac- 
cretion front in the absence of diffusion, which we have 
neglected for the reasons given in ^ If unaccounted for, 
this non-physical excess advection into C-rich material 
leads to premature burning of the small amount of hy- 
drogen that precedes the accretion front. Thus, when 
calculating the burning rate and compositional changes, 
the core C/0 is treated as ^'^Mg, which is essentially 
inert in the nuclear reaction network at these tempera- 
tures. For all other calculations, the core remains half 
and half ^^O by mass. 

The code is evolved until any local thermal gradient 
is steeper than the local adiabatic gradient and convec- 
tion sets in, at which point we consider ignition to have 
occurred. We define the ignition mass to be the total 
accreted mass at the onset of convection, A/jgn = A/tjgn, 
taking into account the time required for the accretion 
front to reach the top of our grid prior to the beginning 

^ http://mesa.sourceforge.net/ 

''' As noted in ^ systems do exist where the mass fractions of 
accreted ^^C and^'^N are non-solar while the other metals still 
have their solar values. However, for the sake of consistency and 
convenience, we scale all the metals by the same value in each 
model. 




P[dyne cm"^] 



Fig. 2. — Time series of envelope profiles for a 1.0 Mq WD 
accreting at M = 10"* Mq yr'^ with Z = 5.0 Zq and X3 = lO""*. 
Each profile in each panel is separated by 317 yr. The thick line is 
the final profile just prior to convection, 1988 yr after the onset of 
accretion. Thus, the first profiles shown represent the layer after 
1037 yr of accretion. From top to bottom, the panels show the 
temperature, energy generation rate (in cgs units of erg g~^ 
accreted ^^C and ^''N mass fractions, and ^He mass fraction. The 
solid lines in the third panel show X12 , and the dashed lines show 
X14. The vertical dotted line shows the pressure of the envelope 
base just prior to convection. 




IQia 1019 
P[dyne cm~^] 



Fig. 3. — Same as Fig. |2] except with Z = 0.5 Zq, time between 
profiles of 634 yr, and tign = 6337 yr, so the earliest profiles show 
the layer after 4434 yr of accretion. 

of the code run. 

5.2. Representative results 
Figures [2] and [3] show envelope profiles for a 1.0 M© 
WD accreting at M = 10"* Mq yr"! with X3 = IQ-^. 
Figure [2] has Z = 5.0 Zq, and Figure [3] has Z ^ 0.5 Zq. 
Each panel shows profiles at 4 different times, each sep- 
arated by 317 yr in Figure [5] and by 634 yr in Figure [31 
The thick solid lines in both figures show the envelope 
profile just prior to convection. 
Figure [2] shows a typical p+^'^C ignition. The cause of 
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Fig. 4. — Same as Fig. [2 except with Z = Zq, X3 = 2 X 10"^, 
and tign = 1082 yr, and thus the earliest time shown is 448 yr after 
accretion begins. 

the ignition is the onset of p+^^C burning, which releases 
more heat than can be radiatively transported away. At 
ignition, the accreted carbon mass fraction {solid line in 
the third panel) at the envelope base {vertical dotted line) 
is essentially the same as in the rest of the accreted layer 
because it has not had a chance to deplete before ignition 
conditions are met. Thus, the ^^N mass fraction {dashed 
line in the third panel) in the layer is also unchanged, ex- 
cept near the base, where the small amount of C-burning 
has slightly raised the ^''N mass fraction. 

Contrast this sequence of events with the 0.5 Zq case 
shown in Figure [31 Here, carbon has already been de- 
pleted prior to ignition. At the time of ignition, the car- 
bon mass fraction at the base is several orders of magni- 
tude lower than in the rest of the layer {solid line in the 
third panel), and the difference has been added to the ^^N 
mass fraction {dashed line in the third panel). The en- 
ergy generation profile {second panel) clearly shows two 
peaks in the envelope: the shallower occurs where ^^C is 
burned and depleted, and the deeper peak is due to full 
CNO cycle burning. It is at the deeper location that enuc 
begins to spike and convection occurs. Since the layer 
must become hotter and denser to burn ^"^N, more mass 
accumulates prior to ignition. The ignition mass in this 
case is 6.3 x 10~^ Mq, 3 times higher than the C-triggered 
case, which has Mjgn = 2.0 x 10~^ M©. 

The disparity in Mign is even greater when compared 
to a case with a non-negligible amount of '^He, because as 
shown in f|4l ^He always begins burning prior to proton 
captures onto ^^C. Figure[3]shows the outcome of a model 
with X3 — 2 X 10^"^. This ignition is triggered by "^He- 
burning after only 1.1 x 10^^ Mq has been accreted. 

5.3. Ignition masses 
In this section, we show the resulting Mjgn for a 
range of M, M, and composition. Figure [5] shows 
Mign vs. accreted metallicity with X3 = 10^^, and 
Figure [6] shows Mign vs. "^He mass fraction with so- 
lar metallicity. The dotted, dashed, and solid lines 
represent M — 0.6, 1.2, and 1.35 Mq WDs, respec- 




FlG. 5. — Ignition mass vs. metallicity for X3 = 10 * and dif- 
ferent WD masses and M's. Dotted, short-dashed, and solid lines 
are for M = 0.6, 1.2, and 1.35 Mq WDs, respectively. The top 
(higher Mign) curve for each WD mass has M = 10~^ Mq yr~^, 
and the bottom (lower Mign) curve has M = 10~* Mq yr"-*^ for 
M = 0.6 Mq and M = lO"'^ Mq yr''^ for M = 1.2 and 1.35 Mq. 
The long-dashe d lin e shows the analytic scaling of Mign oc Z~^^^, 
as discussed in i|5.3l 



10-3 




Fig. 



Same as Fig. [5] except vs. X3. The metallicity is the 



solar value. The analytic scaling from fi5.3\ of Mign 
displayed. 



(X X. 



-2/5 



lively. The top (higher Mign) curve for each WD mass 
has M = IQ-^ Mq yr'^ The bottom (lower Mign) 
curve has M ^ 10"^ M© yr^^ for M = 0.6 M© and 
M = lO"'^ M© yr-i for M = 1.2 and 1.35 Mq. The 
highest M shown for M = 0.6 Mq differs from that of 
the other masses because a 0.6 M© WD accreting at 
10~^ M© yr~^ burns hydrogen stably without a ny no- 
vae dS ienkiewicd 119751 fl980: Paczvnski 1983; Fui imotd 
ll982bHShen fc Bildsten 20071 : INomoto et al.l l2007) . 

Some general trends are clear. Increasing the WD mass 
and accretion rate lower the ignition mass. A larger mass 
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is equivalent to a smaller radius, both of which contribute 
to a higher value of g, resulting in a higher temperature 
and density for a given envelope mass, and thus ignition 
conditions are reached at lower envelope masses. Higher 
accretion rates translate into higher compressional lu- 
minosities and higher temperatures, also resulting in ig- 
nition for smaller envelope masses. The basic trend of 
lower ignition mass with higher values of ^^C and ^He is 
also sensible: more fuel means more burning and quicker 
buildup to ignition. 

However, as we have discussed previously, changing the 
composition does not just lead to a change in the energy 
generation rate. If the ^He or ^^C is depleted, conditions 
for burning the next fuel will have to be reached for igni- 
tion. This is the cause of the inflections that are present 
in most of the Mjgn-curves in FigureO and for the abrupt 
changes in Mig„ that are seen in Figure [HI Once a crit- 
ical ^^C- or '^He-mass fraction is reached, the ignition 
becomes qualitatively different. Note that carbon deple- 
tion does not have as drastic an effect as '^He depletion 
because depleted carbon is converted to another burning 
catalyst, ^^N, whereas ^He is depleted to hydrogen and 
^He, negligibly increasing the amount of hydrogen fuel. 

The long-dashed lines in Figures [5] and [6] show the 
analytic power-law scaling of Mjgn with Z and as 
derived from equation ([7]). Assuming that the con- 
stant pressure logarithmic temperature dependences of 
the ^^C- and ^He-burning rates are Xi2 — X3 — 15 yields 
Afign oc Z"^/^ for ^^C-triggered novae and Mign oc X^'^^^ 
for "^He-triggered novae. The difference in the exponent 
is due to ^He burning via self-reactions, as opposed to 
^^C burning via proton captures. The numerical results 
of Figure [5] match the analytic scaling fairly well. The 
high X^ models of Figure [6] also match the analytic ex- 
pectation well. However, for lower values of X3, ^He is 
depleted and does not release enough heat to trigger the 
nova. These novae are instead triggered by ^^C or ^''N, 
and since the metallicity is constant along each curve, 
Afjgii is also roughly constant below a critical value of 

Figure [7| shows lines of constant Afign as a function of 
M and M {solid lines) for Z ^ Zq and X3 = 10^^, 
obtained by quadratically fitting the numerical results 
detailed above. The contours are evenly spaced with 
0.2 dex of Afign between each line. The thick solid lines 
have Afign = 10~^ and 10"'' Mq as labeled. Also shown 
are contours of constant recurrence time {dotted lines), 
labeled in years. The hatched region w here hydrogen 
burns stably in a steady state is taken from lNomoto et ahl 
(|2007( ). It is clear that a recurrent nova system accret- 
ing solar metallicity material with low "^He mass fraction 
must have a WD > 1.2 Mq to achieve a recurrence time 
< 30 yr. 

Our ignition masses are a facto r of 1.5 — 3 times h igher 
than the analytic estimates of iFuiimotol (ll982bD an d 
1.1-2 times lower than the work of 'Ma cDonaldl (| 19831 ). 
Given the differences in method and updates in equations 
of state and opaci ty, these disc repancies are not signifi- 
cant. The work of lYaron et al.l (|2005i) assumes diffusive 
accretion onto initially bare C/0 cores, which invalidates 
our assumption of accretion onto C-poor material, so the 
fact that their Afjgn's are 2 — 4 times lower than ours is 
not surprising. A comparison of our X3 = 10^^ results 
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Fig. 7. — Lines of constant A/ign in M — M space for Z = Zq 
and X3 = 10~* (solid lines). Contours are evenly spaced with 0.2 
decades of Mign between each Une. Thick solid lines have Mign = 
10~^ and 10~* Mq as labeled. Also shown are lines of constant 
recurrence time [dotted lines), labeled in yea rs. The range of stabl e 
H-burning (hatched region) is obtained from lNomoto et al.l 1)20071 '). 

with those of |To wnslcv fc Bildsten' (2004) (and Fig. 1 of 
iTownslev fc Bildsten 2005 ) reveals a large discrepancy, 
with our Afign's higher by a factor of 3 — 10. The dif- 
ference is predominantly due to their lack of abundance 
evolution, which eliminates the possibility of "^He deple- 
tion. Theirs is a valid assumption for higher mass frac- 
tions X3 > 2 X 10~^, but is not correct for X3 = 10^^, 
and leads to premature nova ignitions. 

Throughout this study, we have assumed that Afign and 
iign are the accreted mass and accretion time prior to the 
onset of convection. However, the convective zone takes 
some time to grow and reach the surface, at which point 
the nova outburst becomes observable. It is conceivable 
that a significant amount of mass is accreted during the 
convective phase, which would cause we have cal- 

culated to be lower than the actual time between novae. 
To quantify this error, we consider the growth time of 
the convective zone, which is ^growth cpTb/e, where 
the temperature and burning rate are evaluated at the 
onset of convection. For each model, we find that tgrowth 
is always 20—100 times shorter than the time to the onset 
of convection, so that the ignition times and masses we 
have calculated are only smaller than the actual values 
by a few percent, and our assumption is justified. 

5.4. Ignition pressures and the misuse of "Pait" 

IFuiimotol ()1982aD is often cited to support the claim of 
a critical ignition pressure for nova ignition, Pcrit, which, 
depending on the study making the claim, has a value 
between 10^^ — 10^° dyne cm~^. However, the ignition 
pressure is clearly not constant, as we show in Figures 
[5] and [3 These figures are identical to Figures [5] and [H] 
but with Pign along the ordinate axis instead of Afign; the 
ignition pressure is related to Afign through equation 
The ignition pressure varies by a factor of 100 from one 
extreme of high accretion rate, mass, and ^"3 to the other 
extreme and is not constant. Moreover, the or iginal pa- 
per actually makes no such assertion. Instead, IFuiimotol 
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a range of M and M. Dotted, dashed, 
M = 0.6, 1.2, and 1.35 Mq, respectively, 
each WD mass have M = 10~^ Mq yr~^. 
for M = 0.6 Mq has M = 10"** Mq yr" 
curves for AI = 1.2 and 1.35 Mq have 10" 
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Fig. 9. — Same as Fig. \E\ except vs. X3 with solar metallicity. 



(jl982af l states that there is a minimum ignition pressure 
necessary to produce a strong nova-like outburst powered 
by hydrostatic shell expansion. The ignition p ressure is 
left as a free parameter and is calculated in iFuiimotol 
(fl982b.l as a function of WD masses and accretion rates. 

5.5. Pre-ignition luminosities 

The depletion of fuel wiU increase the WD surface lu- 
minosity above the compressional luminosity. Such an 
increase in the pre-nova light curve would be a sign that 
the carbon abundance is not enhanced, or that the ^He 
abundance is non-negligible. There are several energy 
scales to keep in mind: the accretion, thermal, ^^C- 
burning, and '^He-burning energies per mass are given, 
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Fig. 10. — Surface luminosity, excluding the accretion luminos- 
ity, and effective temperature as a function of time for a 1.0 Mq 
WD with M = 10~* Mq yr"i, X3 = 10"*, and various metallic- 
ities. The metalhcities shown are 0.1, 0.2, 0.5, 1.0, 2.0, & 5.0 Zq, 
increasing from right to left. 



respectively, as 



therm 



GM 

kT 
firUp 



■==1.3 X 10 



17 Ml erg 
R9 g 



a.4x IQis 



^12^^12 =8.8 X 10 
X3E3 = 2.1 X 10 



14 



X12 erg 



10-3 g 

15 ^3 erg 
10-3 g ' 



(11) 



The luminosities associated with these energies can be 
obtained by multiplying them by the accretion rate, 
with an additional factor / — 1.75 for the ther- 
mal/compressional luminosity. The radius is typically 
a few X 10* cm, so the accretion energy is roughly two 
orders of magnitude larger than the thermal or nuclear 
energy available during the accretion phase. However, 
the accretion luminosity from a disk is variable and 
greatly reduced in disk quiescence. Thus, it is still 
possible to observe the luminosity produced from the 
interior of the accre t ed en velope, as demonstrated by 
iTownslev fc BildstenI (|2003f ) in their work on relating ef- 
fective temperature, Tes, to M. 

Figure [10] shows the surface luminosity, ignoring the 
accretion luminosity, and effective temperature as a func- 
tion of time for a 1.0 Mq WD with M = lO"* M© yr'^ 
X3 = 10-'', and metalhcities of 0.1, 0.2, 0.5, 1.0, 2.0, & 
5.0 Zq, increasing from right to left. The 0.5 Zq model 
is interesting in that there is enough carbon to power 
significant nuclear luminosity, but not enough to trigger 
a nova. For this case, the compressional luminosity is 
< 0.5 Lq and T^s < 5.5 x 10* K for - 50% of the nova 
cycle. Four thousand years after the onset of accretion, 
carbon is burned and depleted, causing the luminosity 
from inside the envelope to double within a span of only 
1000 yr. The burning of ^''N begins - 1000 yr after C- 
depletion, and the luminosity and effective temperature 
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Fig. 11. — Same as Fig. 1101 but with Z = Zq and varying 
X'i. The ^He mass fractions shown are 0.1, 0.2, 0.5, 1.0, 2.0, & 
5.0 X 10""^, increasing from right to left. 

rise to 1.3 Lq and 8 x 10"' K just prior to the CNO cycle- 
triggered ignition. 

Changing also has an effect on the surface lumi- 
nosity and Toff. Figure [TT] shows the surface luminosity 
and effective temperature as in Figure flOl but with fixed 
metallicity Z = Zq and varying ^He mass fractions of 
0.1, 0.2, 0.5, 1.0, 2.0, & 5.0 x lO'^, increasing from right 
to left. Here, the interesting case is X3 = 10~^, which has 
enough '^He to produce significant energy, but not enough 
to trigger the nova. The surface luminosity reaches Lq 
and Teff rises to 6.5 x 10'' K after only 1500 yr of accre- 
tion, but the surface then cools as the envelope succeeds 
in depleting the accumulated ^He. The surface bright- 
ens when C-burning commences, and again dims when 
carbon is depleted, finally rising to 2 Lq and 8 x lO'' K 
when ^^N begins burning. 

6. CONCLUSIONS 

Motivated by uncertainties in classical nova core- 
mixing mechanisms and the lack of metal enhancements 
in some nova ejecta, we have quantified the effects of 
composition on nova ignitions (see Fig. O under the 
assumption that the underlying material is C-poor and 
diffusion thus unimportant, as appropriate for accretors 
with large helium buffers or 0/Ne cores. We have found 
that for carbon mass fractions < 2 x 10^"^, ^^C is depleted 
and converted to ^''N without releasing enough heat to 
trigger a nuclear instability. The layer continues to ac- 
crete until '^^N can capture protons, leading to a nova 
triggered by the full CNO cycle and an ignition mass 
larger than the carbon-ignited case. Thus, the ignition 
mass increases by a factor of ~ 3 as the metallicity is 
decreased from 5.0 Zq to 0.1 Zq. The critical carbon 
mass fraction is near-solar and is thus relevant to sub- 
solar metallicity systems as well as systems with evolved 
secondaries that have undergone CNO processing of ^^C 
to ^"^N. We have also examined the effect of accreted ^He 
(see Fig. |6]), which can reach mass fractions of 4 x 10"'^ as 
an evolved donor's interior is uncovered by mass transfer. 
For X3 > 2 X lO^'^, '^He-l-^He reactions trigger novae with 
Affgii a factor of ~ 3 times smaller than the C-triggered 



case. 

The dependence of Mign on accreted composition will 
affect population-averaged nova rates: naively, high- 
metallicity environments would have nova rates higher by 
a factor of ~ 3 than systems with sub-s olar metallicities 
(such as novae in globular clusters; iShafter fc Quimbvl 
I2OO7D or evolved donors that have undergone CNO pro- 
cessing. However, the existence of '^He would have to 
be taken into account, because old systems with donors 
that have undergone significant mass loss could have 
X3 > 2 X 10~^ and would thus have high nova rates, re- 
gardless of the accreted metallicity. A proper prediction 
of the effect of donor composition requires a population- 
synthesis calculation that includes further complications 
such as binary and donor evolution. We leave this exer- 
cise for future work. 

To date, most observations only report galactically- 
avera ged nova rate s, altho ugh some M31 stud- 
ies diCiardullo et all 119871 : ICapaccioH et all Il989t 
iDarnlev et al.l [2006( 1 have found that M31's bulge pro- 
duces more novae per stellar luminosity than its disk by a 
factor of ~ 5. On the other hand, galactically-averaged 
nova rates do not see any morphology dependence, 
finding instead that the luminosity-specific nova rate 
(LSNR) is roughly constant across all galaxy types at a 
value of 2 ± 1 yr^^ (10^° Lq,k)~^ (Williams & ShaftiB 
120041 ) ■ where Lq^k is the K-h&wd solar luminosity. The 
LMC a nd SMC (and possibly Virgo dwarf elliptical 
galaxies: [Neill fc Sharal [2005l. which have LSNRs higher 
by a factor of 3, are exceptions. These measurements 
have large error bars due to small number statistics 
and issues of completion caused by both extinction and 
infrequent observations. A better measurement of nova 
rates will come with new deep optical surveys with 
high cadences such as Pan-STARRS-1, Pan-STARRS-4, 
and the Large Synoptic Survey Telescope, which will 
see thousands of novae every year. These will reduce 
the nova rate error bars and also possibly allow us 
to measure rates in different populations within other 
galaxies besides M31. 

In addition to this population-averaged observable, the 
composition could also have a detectable effect on indi- 
vidual systems. In particular, the depletion of fuel can 
significantly increase the surface luminosity above the 
baseline set by the entropy released during compression 
of the accreted layer (see Figs. [TOl and [TT|) . While this 
increase is still well below the accretion luminosity as- 
sociated with gravitational energy release, it would be 
visible while the system was in disk quiescence. Recur- 
rent novae, in particular, would be ideal systems in which 
to observe this brightening in quiescent luminosity due 
to their short recurrence times of '~ 30 yr. 

These observables are dependent on the assumption 
that convection is not initiated in C-enhanced material. 
If instead H-rich envelope material penetrates into C-rich 
material and triggers convection there, the accreted com- 
position will have little effect. Thus, if these composition- 
dependent effects are observed, they will provide evi- 
dence that convection for many novae is initiated in C- 
poor material, and that CNO enrichment for these novae 
is due to convective shear mixing or overshoot. 
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